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CL,; Abstract 

r-| ! The chiral expansion of the tttt amphtude to the order of two loops was expressed in terms 

of six independent parameters in a previous paper: four of these are shown here to satisfy 
• 1-^ . 

^ ■ sum rules. Their derivation, where crossing symmetry plays a key role, is explained. Their 
d • convergence properties are studied and their practical evaluation, in terms of the available data 
on TTTT phase shifts above 0.5 GeV, is discussed. Below 0.5 GeV, the chiral amplitude itself 
is employed, such that the parameters are determined in a self-consistent way. Some care is 
devoted to the estimate of the errors. 
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1. Introduction 

In a previous paper Q ( referred to as I in the sequel), we showed that the vrvr elastic 
scattering amplitude to the order of two loops in the chiral expansion takes the form of an 
analytic expression that depends on six parameters (a, (3, Ai, A2, A3, A4) which are not 
determined by chiral symmetry. In the present work, the possibility of determining these 
parameters from experiment will be discussed. It will be shown that four of them, Ai,...,A4, 
can be inferred with rather good accuracy from existing tttt scattering data at energies ^/s > 500 
MeV. This determination is merely based on very general properties of the scattering amplitude 
and it is completely model independent. Before entering into the details of the subject a few 
comments on the low-energy constants characterizing higher orders in the low-energy expansion 
may be useful. 

Upon extending the program of the chiral perturbation theory (xPT) beyond one-loop 
order [Q, one faces the problem of the proliferation of low-energy constants: at order 0{p^) 
there are more than a hundred new unknown parameters At one- loop level, the predictive 
power of xPT resides in the fact that the same (combinations of) low-energy constants often 
appear in different observables. This fact, of course, reflects chiral Ward identities and unitarity 
and it is encoded in the effective lagrangian. However, at two-loop order, eliminating unknown 
parameters in the same way by comparing different observables becomes problematic in practice, 
though not impossible in principle. In spite of this difficulty, successful two-loop calculations 
do exist P], in which the influence of the new 0{p^) constants remains rather limited, at least 
in a particular kinematical region. 

In general, it seems plausible that at 0{p^) accuracy, each experimentally relevant process 
will be described by its own, specific set of low-energy constants, which have to be determined 
from considerations operating beyond the strict framework of the low-energy expansion. Models 
aiming at such a determination have been proposed, and already applied in the past both at 
the O(p^) P, and 0{p^) ^ levels. They incorporate low-lying meson resonances into the 
effective lagrangian, thereby extending its domain of applicability. The unknown low-energy 
constants are then obtained by fitting the resonance parameters and couplings first, and then 
integrating out the resonance degrees of freedom (expanding in inverse powers of their masses) . 
While this procedure is essentially unique and presumably rather accurate in the case of vector 
or axial- vector resonances 0, it is less convincing as far as scalar exchanges are concerned: i) 
first, incorporating spin resonances into Cejf in a way consistent with chiral symmetry is not 
a unique procedure. For obvious kinematical reasons, the criterion of matching with QCD at 
short distances is less stringent than in the case of spin-1 exchanges. Too many terms survive 
this criterion, which even requires the presence of scalar couplings with a non-minimal number 
of derivatives (which were not considered in P]). ii) Next, it has been suggested that the 
spectrum and couplings of the scalar mesons are subject to particularly strong distortions due 

^ A brief account of the results of the present paper has already been presented in I (sec. 4) . 
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to unitarity and final-state interactions ||^. It may be misleading to describe the 0+"*" channels 
by a set of quasi-real poles or, equivalently, to identify the resonance masses and couplings in 
the extended effective lagrangian with the corresponding entries in the PDG tables [|I^ . For all 
these reasons, a more realistic and less model dependent method of determination of low-energy 
constants is needed, especially in the symmetry breaking sector which is particularly sensitive 
to scalar exchanges. 

In this paper, such an unambiguous method is worked out in the case of elastic tttt scattering. 
It may be viewed as a generalization of the method of chiral sum rules for various two-point 
functions |[Tl|, |12|, which has been already used to determine some 0{p^) constants in connection 



with the reaction 77 vr^vr'^ [|T2| . The method makes use of analyticity, unitarity and crossing 
symmetry to relate the existing experimental information on the vrvr scattering amplitude at 
medium energies, 0.5 GeV < ^/s < 1.9 GeV, to its behaviour in the threshold region, where xPT 
applies. The comparison yields six sum rules for the six unknown parameters a, (3, Ai, A4 
which define the perturbative tttt amplitude up to two loops. Four of these sum rules converge 
rapidly enough to be practically insensitive to the poorly known high-energy part of the tttt 
amplitude. In the absence of information on this part, the remaining two sum rules are of little 
use in practice. For this reason, our method merely serves to fix the constants Ai, A4, leaving 
the parameters a and (3 undetermined. This, of course, is reminiscent of the derivation of the 



Roy equations |T^, for the vrvr scattering amplitude: here, as well, after imposing crossing 



symmetry and fixed-t analyticity, one remains with two arbitrary subtraction constants [1^, |I4 
The parallel with the Roy equations will be further developed in sec. 2, where the derivation 
of the sum rules is presented in details. 

The TTTT amplitude is a rare example of a low-energy observable which is sensitive to the 
mechanism of spontaneous chiral symmetry breaking, in particular to the strength of quark 



antiquark condensation in the QCD vacuum [15, 16|. Indeed, elastic vrvr scattering provides the 



best way of testing the standard postulate, according to which the quark condensate 

5o = - lim <qq> (1) 

TV 

is large enough (c.f. Bq ~ 1.5 GeV ) to dominate the response to the perturbation by light 
quark masses [|l7|- This hypothesis is at the basis of the standard xPT H, which explicitly 
assumes that the GOR relation 2rhBQ = is not violated by more than a few percent 0|. This 
assumption has not yet been tested experimentally, and it is not a necessary prerequisite for a 
consistent chiral expansion: indeed, it was shown in [|l^, |16[ |12| that a systematic, more general 



expansion (G^PT) can be developed, which neither requires Bq to be large, nor the expansion 
of to be dominated by 2mBo (see for a recent review). The two-loop expression of the 



TTTT amplitude displayed in I holds in full generality, independently of any particular mechanism 
of chiral symmetry breaking. The strength of quark condensation merely affects the values of 
the parameters a, (3, Ai, A4, but not the analytic form of the amplitude. The standard and 
the generalized xPT differ in the expansion of these parameters in terms of quark masses, chiral 
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logarithms and the renormahzed (quark mass independent) constants of Cejf- This expansion 
was worked out very recently [|l^ up to chiral order six in the standard xPT. This kind of 
information is particularly useful for the two parameters a and /?, the expansion of which start 
at the lowest Oij?^ order of the xPT, and higher order contributions are expected to remain 
relatively small. In the standard xPT, both a and /5 are predicted to be close to 1 (within 
10 — 20%, say), while in the alternative situation of a small condensate the value of a would be 
2 — 3 times larger, depending on the actual value of the quark mass ratio ms/rh-, which remains 
a free parameter in the GxPT. In all cases the deviation of (3 from 1 should remain small: for 
more details, see I. In contrast, the expansions of the parameters Ai, A2 and A3, A4 start at 
order 0{j)^) and O(p^), respectively. Consequently, these parameters are more sensitive to the 
unknown 0{p^) constants of £e// and their values are thus expected to be more difficult to 
estimate on the basis of a xPT calculation alone. It is gratifying that the experimental vrvr data 
at medium energies allow to determine Ai, A4, and it is significant that the same data leave 
a and (5 undetermined. Our analysis further suggests that the values of Ai, A4 are practically 
insensitive to the strength of quark condensation. We checked that by allowing a to vary in a 
rather wide range (from 1 to 4), the calculated values of the Aj's were only affected at the level 
of a few percent. In particular, the values we obtain should be used within the framework of 
the standard xPT (see sec. 3.5). 



In order to disentangle the large and small condensate alternatives, it would be ideal to be 
able to exploit the variation of the parameters a, /?, Ai, A4 as functions of quark masses and 
chiral logarithms as predicted in the standard and generalized xPT, respectively. Unfortunately, 
this kind of theoretical information is not suitable for a direct experimental test, since Nature 
has already made its choice of quark masses. On the other hand, the forthcoming high-precision 



data near threshold [^, ^ should allow to measure the parameter a, for which the alternatives 
in question differ by a factor 2 — 3. Obviously, determining all six parameters of the two- loop 
amplitude from a fit to the low-energy data appears practically hopeless. For this reason, the 
model independent determination of the parameters Ai, A4 reported in this paper becomes 
of a fundamental importance. Special care is devoted to the discussion of the experimental 
inputs and of the uncertainties in the results (sec. 3). A reliable estimate of the errors is indeed 
crucial to correctly assess the size of the expected theoretical uncertainty in future experimental 
determinations of the critical parameter a. 

2. Derivation of the sum rules 



We start by giving a brief description of the vrvr amplitude to two loops, in order to set up 
the definition of the parameters for which we intend to establish the sum rules. We will neither 
repeat the derivation nor even give the complete formulas, all these details can be found in I. 
The T-matrix element for the process -k""k^ ti'^tt'^ is expressed in a standard way in terms of 
a single function A(s|t, u), symmetric in t, u, and where the Mandelstam variables s, t, u satisfy 
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s + t + M = 4M2, as 

T"^'"'^(s, t, u) = A{s\t, + A{t\s, u)5'''5'"^ + A{u\t, s)5'"^5'"'. 



(2) 



At two-loop order, the function A{s\t,u) may be written as a sum of two terms 

^2-'°"^ u) = AP^^slt, u) + u) , (3) 



where the former is a polynomial in the Mandelstam variables and the latter collects the cuts 
of the amplitude. The polynomial is of the third order and satisfies crossing symmetry. It may 
be parametrized as 



AP°\s\t, 



u 



a Ml 
F2 3 



V 3 ^ „ 
+^As - 2Mlf + ^ [(t - 2Mlf + {u- 2MIY 

TT 7r 

+je{s - 2Mlf + ^,[{t- 2Mlf + {u- 2Mlf 



(4) 



which displays the six parameters a, /5, Ai, A2, A2, A4 of the two-loop amplitude. The part 
containing the cuts in the complex plane can be expressed in terms of three functions of a single 
variable Wa{z), a = 0,1,2, which are i) analytic in the whole complex plane in the variable 
z except for a right-hand cut on the real axis [4M^, 00], and ii) asymptotically bounded when 
z 00, such that Wa/z^ — > for a = 0, 2 and Wi/z^ 0. In terms of these functions, 
^cut/^u \ written as 



1 



A'^'Hslt^u) 



327T 



+ 



3 
1 

2 
1 

+ 2 



Wo(s) - W2(s] 



3{s-u)Wiit) + W2{t) 
3{s-t)Wi{u) + W2{u) 



(5) 



The fact that the amplitude to two loops has the general structure described in the above 
equations (|l),(i),(iD was first proved in |jl6l. In the subsequent paper I, the explicit form of 
the three functions Wa{z) was obtained f\ 



Wa{z) = Y.<i^)Kn{z) 



n=0 



where 



Ko{z) 



167r2 



In 1 



1 + Jl-AM^/z^ 



(6) 



(7) 



The functions Wa{z) introduced in I differ from the functions Wa{z) given in by polynomials which are 
defined up to an ambiguity described in appendix B of [Rw. 
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is the Chew-Mandelstam function which is already present at one-loop order, and the other 
four functions Ki, K4 are simple combinations of Kq , (A'o)^ and (-^o)^- The set of functions 
w2{z) are third-degree polynomials in z which are tabulated in I. These polynomials depend on 
four constants, a, (3, Ai, A2 which are the same (up to terms of order O(p^) in the amplitude) 
as those which appear in The dependence is cubic in a and /? and linear in Ai and A2. 
Finally, we note that there is a simple relation among the discontinuities of the functions Wa 
and the discontinuities of the S and P partial waves along the right-hand cut, s > 4M^: 

lmfS{s) = lmWa{s), a = 0,2, lmfl{s) = {s-AM'JlmW,{s). (8) 

This relation is actually used in the construction of the functions Wa in a recursive way, starting 
from the expression of the partial waves at order 0{p'^) and using unitarity to generate the 



imaginary part at a higher order. This technique was used for the first time in ||22[ to obtain 
the one-loop vrvr amplitude in the case M-j^ = 0. 

2.1 Roy dispersion relations 

Due to the Froissart bound, the vrvr scattering amplitude obeys fixed t, twice-subtracted 



dispersion relations which involve three subtraction functions. It was shown by Roy [T^ and 



by Basdevant et al. |T3 that these three functions are entirely determined, up to two constants 



(which, for instance, may be taken as the two scattering lengths Oq and Oq), once crossing 
symmetry is imposed. The key step in the derivation of the sum rules is to equate, in the low 
energy region, the perturbative expansion of the amplitude (||) (^) (||) with the Roy dispersive 
representation. We start by rederiving this representation in a form which will be convenient 
for our purpose. 

Consider the s-channel isospin J = 0, 1,2 amplitudes F^{s,t), defined as 

FO(s,t) = :^{3A{s\t,u) + A{t\s,u) + A{u\s,t)}, 
ozvr 

F\s,t) = ^ {Ait\s,u)-Aiu\s,t)} , (9) 
ozvr 

.2, ^ 1 



F'is,t) = — {A{t\s,u) + A{u\s,t)} 



and form a 3- vector 



Fis,t) 



( F%s,t) \ 
F\s,t) 



(10) 



The Froissart bound allows one to write a fixed t dispersion relation with two subtractions 

F(s,t) = C.4a+(t) + (.-«)a_(t)] + - / — + ^ ImF(x,t) , (11) 

vr J AMI \ x — S X — u I 
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where Cst, C'su and Cm are the s — t, s — u and t — u crossing matrices 





( 1/3 


1 


5/3 \ 




/ 1/3 


-1 


5/3 ^ 




(I 





o\ 


Cst — 


1/3 


1/2 


-5/6 




-1/3 


1/2 


5/6 


) Ctu = 





-1 







I 1/3 


-1/2 


1/6 j 




I 1/3 


1/2 


1/6 ) 









1 / 



• (12) 



We have also introduced the vectors 



a+(t) 



/ ao(t) \ 


V «2(t) ; 



a_(t) 



( ° 1 

ai(t) 

V ; 



(13) 



which collect three arbitrary functions of t. As was observed in [III, M , the constraint of crossing 
symmetry reduces this arbitrariness down to two real constants. From (|Tl|), F manifestly 
satisfies s — u crossing by construction, i.e. F{s,t) = CsuF{u,t). One must impose s — t 
crossing symmetry, and then t ~ u crossing symmetry will follow automatically. 

It proves convenient to split the integration range into two parts, [4M^, E^] and [E'^, oo]. In 



the classic study of the Roy equations |23[ E was taken to be rather large, E ~ 1.5 GeV. Here, 
for the purpose of generating a form which can easily be compared with the two-loop expression, 
we will take E to be sufficiently small, so that in the range 4M^ < x < E^ the contributions 
of the partial waves with £ > 2 to the imaginary part of the amplitude are negligibly small as 
compared to the S- and P-wave contributions. This condition is satisfied for E 1 GeV. For 
X < E"^ we can write 



lmF{x,t) 



( Im/o°(x) \ 


V Im/2(x) J 



3(1 + 



2t 



X 



AMI 



I \ 

Im/i(a;) 




(14) 



Inserting this into the dispersion relation ([TT|), there will appear three functions W^{z) which 
are analogous to the functions Wa{z) in the two-loop expression: 



W^{z) = — dxhnJSix) 

° TT JiMi x'^ X — Z 



(a = 0,2) , 



Wf{z) 



We collect these functions into vectors 



Wliz 



( <W \ 


V wiiz) ) 



W^iz) 



dx \mfl{x) 
TT Jam^ x{x — 4:M^) X — z 



( ° ^ 

W^{z) 

V ) 



(15) 



(16) 



The amplitude F(s, t) can then be expressed as a sum of two terms 

F{s,t) =W^{s,t) + Z^{s,t) . 



(17) 
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The first term is constructed from the integrals in order to satisfy crossing symmetry 
exactly: 

+Cst[W^{t) + 3(s - u)W^{t)] . (18) 

The second term, Z^{s,t), collects what is left from the original dispersive representation. 
Remarkably, it can be expressed solely in terms of integrals over the high-energy range [E"^, oo] 
and of three arbitrary functions of t. This is achieved by redefining the original subtraction 
functions a_|_(i(:) and a_(t) into new ones h^{t) and b_(t) (which are equally arbitrary for the 
moment), collecting into them as many terms as possible: 

b+(t) = a+(t) - W^{t) + 3(t - 4M^)(t - 2M2)(1 + QJ- / ^C.Jm W^(x) 

n J AMI 

1 fE^ fir 

MO = a.(0 - m^(t) - S(t - 2Af,)-l^,^ ___,,„;;(,) 

-it-AMD- ^-^^^ImF(^-i)(x,t) . (19) 
n Je^ x^{x - AMI) 

We have introduced the amplitudes with given isospin in the t-channel which are defined, as 
usual, by 

F^'^=^\s,t)=j2{Cst)aiF\s,t) . (20) 

These explicit relations between the old and the new subtraction functions will only be useful 
at a later stage for the derivation of slowly convergent sum rules for the parameters a and jS. 
Z^{s,t) has now the following expression: 

Z^(s, t) = Cst{h+{t) + (s - M)b„(t) + us[h+{t, us) + (s - u)h_{t, ms)]} , (21) 

where the vectors h^{t,us) and h_{t,us) are formed in analogy to a±, W± (see (PTSD, (|TBp) 
from the three functions 

_ -1 .oo dxi2x + t-AMDlmF(^^=^\x,t) 
na[t, US) - ^ j^^ ^i^ + t- AMl)[x^ + x{t - AMI) + us] ' ' ^ ' 

and 

c/x Im F(^*=i) (x, t) 



h^{t,us) jj,, ^i^^t_4^Ml)[x^ + x{t-AMl)+us]' ^^^^ 

As announced, these dispersive integrals extend over the high energy range [i?^, oo\. 

Up to this point, E was bounded only from above by the requirement that the contribution 
from the D-waves can be neglected for energies smaller than E. We can now also take E to be 
sufficiently large compared to twice the pion mass. In other words we require 

AMI «E'^ <1 GeV^ . (24) 
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Then, in the low energy region, the ratios s/E"^, t/E'^-, u/E"^ will be of order 0{p^) in the chiral 
counting; consequently, we can expand the integrals ha(t,us) in powers of these quantities and 
drop the terms which are of chiral order 0{p^) or higher in the amplitude (this is legitimate 
since we intend to equate this form of the amplitude with the expansion to two-loop chiral 
order). Hence, we can replace in Eq. (|2T|) 

ha{t,us) = ha{0,0) + tdtha{t,0)\t=o + 0{p^) (a = 0,2), hi{t,us) = h{0,0) + O{p^) . (25) 

Finally, we must impose s — t crossing symmetry (to the same chiral order) on Z^{s,t), i.e. 

Z^(t,s) = C,iZ^(s,t) + 0(p«) . (26) 



The standard trick of setting s = in ( pG]) shows that the subtraction functions ba{t) must 
be polynomials once the integrals ha{t,us) have been expanded (^Sf). These polynomials are 
entirely determined up to two constants, a and b: 



/\/f2 1 

bo{t) = ^(5a - 86) + 2tb - -t{t - AMI) 



ho + Q{t - 2Ml)hi + 5/i2 



hi{t) 



t r 
6 



2ho + 6{t - 2Ml)hi - 5/i2 



2M^ t r 

— ^(a + 26) ~tb- -{t - ml) 2ho - 6(t - 2Ml)hi + /12 

3 O L 



where we have introduced the notation 



K = ha{o,o), K = dMt,o)\t=o . 



(27) 



(28) 



The crossing relation (^) is then identically satisfied provided the following combination of 
high energy integrals vanishes in the chiral limit: 



(29) 



(Using resonance saturation, this constraint may be converted into an amusing relation between 
the masses and the widths of the /2(1270) and the p3(1690) resonances ^) This completes the 
construction of a dispersive representation of the amplitude F(s,)f:) in terms of two arbitrary 
parameters, a and b. Explicitly, from ([l^) and (|27D one obtains the following expression for the 
corresponding function A{s\t,u) 



^A^^^Pis\t,u) 



+ l(s-u)wfit) + mfit) 



+1(3 - t)wf{u) + iwiiu) 

+hi{-s^ + 6My - 8M^s) + lih'o - h'^)stu - i/i2s(s - AM^) + i(/io - h2)tu 

+bs + lia - 4b)Ml +0(/) . 

(30) 

^ This relation reads: 63TpMj = 5TfMp, where Mf and Mp are the masses of the /2 and the mesons 
and Tf, Tp are the respective tttt partial widths. 
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Expression (^) represents the final form of the dispersive nn amphtude. We recall that the 
functions are computed from the imaginary parts of the partial- wave amplitudes /o(s), 
fl{s) and /o(s) (see (0)) in the domain 4M^ < s < E"^, i.e. in the region where experimental 
information is incomplete, while the constants ha and h'^ (|2^) can be evaluated using existing 
experimental data at higher energy: the part involving these functions corresponds to the so- 
called driving terms in the Roy equations |2^. Indeed, from A'^^'^^{s\t,u) one can compute the 
real parts of /g , fi and /g , which are thus given in terms of the corresponding imaginary parts 
and the driving terms. These integral equations, together with the elastic unitarity relations 
(which hold to good accuracy below the KK threshold) constitute the Roy equations. Numer- 
ical solutions to these equations were constructed in p3| ||2^ taking as the two free parameters 
the scattering lengths af] and Oq. Constraining these solutions to match with the experimental 
data in the 600 — 800 MeV region, it was shown that is left essentially undetermined but 
that Og gets strongly correlated with Og (this relationship is known as the Morgan-Shaw band 
in the literature p^). We are going to show that in a large domain of parameters a and f3, 



the perturbative amplitude A'^~''°"'P{s\t,u) of Eq.(^ provides a rather accurate representation 
of the numerical solutions of the Roy equations (up to energies y/s ~ 0.5 GeV), provided the 
values of the parameters Ai, A2, A3 and A4 are properly chosen. 

2.2 Equating the dispersive and perturbative amplitudes 

Let us now restrict ourselves to the domain where the Mandelstam variables s, t, u are very 
small (compared to 1 GeV^). In this domain, A'^'^^^{s\t, u) must be identical with A^~^°°^{s\t, u), 
except for terms of chiral order 0{p^) or higher. Let us then subtract the part which contains all 
the cuts of the two-loop expression, A'^'^^{s\t, u) (see (|^)) from the dispersive representation (|30|) 
of A{s\t,u): this difference must be identical to the polynomial part of the chiral amplitude, 
AP°\s\t,u) (c.f. Eq.(§). Indeed, at low energies, the discontinuity in the differences 



Im 



W^{z)-Wa{z) =Im/o"(z)-Im/o"(^)|2-,,,, (a = 0,2) (31) 



is, by definition, a quantity of chiral order eight. A similar relation also holds for the isospin 
one case. In the range where z « 1 GeV^, these discontinuities can be neglected and the 
differences W^{z) — Wa{z) are analytic functions which are well-approximated by polynomials 

32vr[W^f (z) - Wa{z)] = i: ly + 0{z') . (32) 

A:=0 

Inserting these expansions into the difference A'^''P{s\t, u) - A^-^°°v{s\t, u), we generate six sum 
rules by recognizing that the polynomial obtained in this way should vanish up to terms of 
order 0{p^). It is then a simple exercise to derive: 

J, = lil'o - ID - 3/i + ^Mlill - II - + ^-fiho - Ah,) 
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3 
2" 



{ho - h2) 



(33) 



As 



327r 



.Xll-ll) + ll + —{h',-h',-9h, 



f! 2^ 2 



9 



where the entries are defined in (|3^). 



The remaining two equations similar to (|33|) concern the parameters a and p. However, 
these depend exphcitly on the two subtraction constants a and b of the Roy representation (pOT): 



a 



JL 



327r 



2Ml{ho + 2/12) + -Mt{K -K- + a 



(34) 



l2Mt { hi - hi - hi ) + 327r 4M>i + h 



Consequently, one cannot, a priori, make use of these equations to determine a and j3. 

The sum rules (|3]) and (|3^) represent the minimal set of necessary and sufficient conditions 
ensuring the compatibility of the dispersive representation ( |30D with the chiral expansion of the 
amplitude up to and including the two-loop accuracy. It is worth stressing the role of crossing 
symmetry in the derivation of these sum rules. 

2.3 Barely converging sum rules for the parameters a and (3 

The convergence of the high-energy integrals in Eqs. (|33|) is guaranteed by the Froissart 
bound. Actually, as will be seen in sec. 3, the sum rules (|33| ) are rapidly convergent and rather 
independent of the details of high-energy asymptotics. Additional information on the high- 
energy behaviour of the amplitude would yield additional sum rules. The standard Regge pole 
phenomenology, for instance, suggests that the /t = 1 t-channel isospin amplitude F^^*"^^(s,t) 
behaves asymptotically as s"''*^*-' where dpit) is the p trajectory with the intercept ap(0) ~ 1/2. 
According to this picture, the amplitude F^^^^^'^ satisfies a once-subtracted dispersion relation. 
This idea was used in the past to generate a sum rule for the P-wave scattering length a\ and 

Comparing the once-subtracted with the twice-subtracted 



for the combination 2aQ 



Sag 



fixed t dispersion relation ([TT|), one obtains the following expression for the subtraction function 
ai{ty. 

(35) 



71 J AMI 



Similarly, the amplitude with F^^' '^\s,t) is dominated by exotic t-channel exchanges and it is 
expected to be asymtotically suppressed. It is even conceivable that it satisfies an unsubtracted 
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dispersion relation ||2^. Under this assumption, one obtains the following sum rule for the 
subtraction function 02 (t) in (|ll|): 

1 f°° dx 

vr J4M2 



1 rco fjrp 

a2{t) = - — {2x + AMI - t)ImF(^*=2)(x,t) . (36) 

vr J4M2 



The integral (|36|), if convergent at all, converges even more slowly than the sum rule (^S]). The 
asymptotic behavior of F^^*"^-* is expected to be dominated by s^^p"^ arising from the p — p 



Regge cut ||2^. Assuming the validity of Eqs. ( p5D and (pGf ) at t = 0, one can use them in ([19|) 



and then, from (^) express the two subtraction constants a and h\ 

b = 7^ r\ {2/o°(^) - + ^jiT^/i (^)} 

GtT J4M2 X"' X — 4M^ J 

+- r ^-^^T7?vImF(^*=^)(x,0), (37) 

(a + 26)M^^= ^ / ^ -Im{(l + ^)[2/°(x) + /o^(x)]-9^^/,^(x)} 
27r J4M2 X X X — 4M^ J 

+ - r ^(x + 2M2)ImF(^*=2)(x,0). (38) 
vr J£:2 X 



These relations promote Eqs. (|34D to two additional sum rules for a and /3. 

The slow convergence of the high-energy integrals in Eqs. (^) and p8| ) prohibits any 
practical use of these sum rules for a and p. The sum rule for /5, which is expected to converge 
better, leads to values /3 = 1.2 — 1.4 which have the correct order of magnitude (as compared 
to the xPT prediction), but the uncertainty due to the high-energy tail of the integral in ( P7| ) 
is difficult to estimate. Needless to say, the evaluation of the sum rule for a which involves the 
integral ( ^Sf ) is even more problematic. 

In the following, these barely convergent sum rules for (3 and a will be ignored. The 
above discussion should mainly serve as an explanation of why, in contrast to Ai,...,A4, the 
two parameters a and f3 cannot be determined from the existing vrvr scattering data at medium 
energies. New, high precision low-energy experiments are, indeed, unavoidable for this purpose. 

3. Evaluation of the sum rules 

In this section, the details of the evaluation of the four sum rules (|33|) are presented. The 
input data above 0.5 GeV and their use are discussed in subsection 3.1. The treatment of 
the low energy part is explained in subsection 3.2, where the final results are also presented. 
Subsection 3.3 is devoted to a careful analysis of error bars. Finally, subsections 3.4 and 3.5 
contain a comparison with related works, in particular with the recently published standard 
XPT two loop calculation [0. 

3.1 Energy region above 0.5 GeV 

The main sources of experimental information concerning the vrvr amplitude are the produc- 
tion experiments (see ^^), which give reliable results in the range 0.5 < i/i < 1.9 GeV. These, 
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together with (p4D, imply that the terms hi and h[, which are integrals over the range [E"^, oo], 
can essentially be evaluated using experimental data. The contribution from the asymptotic 
region y/s > 1.9 GeV can be estimated using Regge phenomenology and it turns out that this 
contribution is fairly small. This is illustrated in table 1 below, which shows the contributions 
of the various energy ranges to these integrals. 



v^GeV 


[0.5,0.95] 


[0.95,1.9] 


[1.9, cx)] 


total 




-6.15 


-0.59 


-0.08 


-6.82 


/i2(GeV-^) 


0.17 


-0.15 


0.00 


0.02 




-4.34 


-0.11 


-0.002 


-4.45 


/i'o(GeV-^) 


8.15 


-0.54 


-0.13 


7.48 


/i'2(GeV-^) 


11.82 


-0.20 


0.00 


11.62 



Table 1: Contributions from the various energy ranges to the integrals hi and h'i, with 

E = 0.5 GeV 

The table shows clearly that the integrals are rapidly convergent and are dominated by the 
region ^/s < 1 GeV with the exception, however, of h2. In that case, the integrand is very 
small below 1 GeV because of cancellations among the three isospin contributions, and the whole 
integral is much suppressed as compared to the others. This is reminiscent of the approximation 
made in [^ of setting the combination of amplitudes with It = 2 equal to zero in sum rules 



derived there for Ai, A2. One should be aware that this simple approximation breaks down for 
/ig; as can be seen from the table, /ig is even larger than h^. 

We now specify in more detail how we have treated the experimental data in order to obtain 
the numbers given above and the central values for the Aj's as given in the sequel. Let us first 
discuss the region below 1.9 GeV. The region was divided into two subregions which were 
treated somewhat differently: 

a) 0.95 < < 1.9 GeV and b) 0.5 < < 0.95 GeV 

• In the subregion a) we have employed the data obtained by Hyams et al. ||2^ (based 
on the CERN-Munich production experiment [Q, which has by far the best statistics 
to date) in the form of the analytic K-matrix type parametrization which the authors 
provide Phase shifts and inelasticities have been determined for partial waves up to 
£ = 3, which means that the contributions from the resonances /2(1270) and p3(1690) are 
automatically included. Partial waves with £ > 4 were ignored: including the resonance 
/4(2040) in the narrow width approximation proved to have a negligibly small effect. In 
the case of the isospin 1 = 2 and the partial wave i = 0, we have used the data of 



Hoogland et al. [^ which have better statistics in this channel, and we have neglected 
the contributions of the higher partial waves. 



There are a few obvious misprints in the pubHshed formulae. 
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In the subregion b), some care is needed because the results are very sensitive to the 
values of the phases there, particularly to the S and P partial waves. Concerning the 
/ = S-wave, several slightly different analyses of the CERN-Munich data have been 
performed. We shall use, as before, the data of Hyams et al. which extend down to 0.6 
GeV, and also the results of Estabrooks and Martin who provide the phases down to 
0.5 GeV. We have performed a fit to both sets of data using the parametrization proposed 



by Schenk |^| 
iaii5l{s) = 



4M2 /AMI 



P 



So 



So 



P 



AMI 



(39) 



where c, d and sq are three parameters to be determined by the fit. In principle, one 
could take a\ to be a free parameter as well, since we intend to use ( P^D rather far from 
the threshold, in the energy region from 0.5 to 1 GeV. However, we have preferred to keep 
with the idea of ||33|, taking to be the experimental value of the scattering length, in 
order not to give too much weight to the first few points around 0.5 GeV. The quality 
of the fit is practically independent of the exact value of a\, within the range allowed 
by experiment. The numerical values that we have used are collected in table 2. In the 
case of the P-wave, the data deduced from the CERN-Munich experiment are known to 



be in slight conflict with the results based on the Roy equations below 0.7 GeV ||3^, |35 
In that case, the Roy equation studies suggest that the 5\ phase is very close to a pure 
Breit-Wigner shape below the resonance. Based on this prejudice, we have again used 
the parametrization (|39D , imposing the experimental value of the scattering length and 
adjusting the remaining three parameters in order to fit the experimental values of the p 
mass and of the p width, and such that = 0.95 GeV) matches to the experimental 

value of |2^. If one uses a real fit to the data of p9| instead, then the difference in 
the results would be of the order of 5%, which means that they are perfectly compatible 
within the errors. Finally, for the I = 2 i = partial wave we have fitted the parameters 
of the representation ( ^9]) to the data of |31[]. In that case, we have used the same 
parametrization in both region a) and region b). 



For comparison, we have also used the parametrizations of the / = S-wave given by |3^, |37 
of the 1 = 1 P-wave given by [^] and the 1 = 2 S-wave given by |29|. Our final results are 
consistent with these fits to the data. 





4 


c 


d 


^(MeV) 


1=0 e = o 


0.26 


0.2535 


-0.0200 


843 


1=2 i = 


-0.028 


-0.2325 


-0.0160 





1=1 £ = 1 


0.038 


0.2560 10-3 


-0.600910"^ 


769 



Table 2: Numerical values of the parameters used in for the S and P partial waves. The 

values of the scattering lengths are taken from 
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In the asymptotic region we have assumed that the amphtudes are given by Regge phe- 
nomenology, and we have used the same parameters as in formulae (17), (18) and (19) of |]40[ . 
This representation was assumed to hold for ^/s > 3 GeV, and in the region between 1.9 and 
3 GeV the amplitude was determined by linear interpolation. 

3.2 Energy region below 0.5 GeV and results 

We must now deal with the terms in the sum rules. According to the definition (|32D , 
they are obtained by making a Taylor expansion of the difference W^{z) — Wa{z). Here, Wa{z) 
is known analytically in terms of a, (3 and is a linear function of Ai, A2. Since E is constrained 
not to exceed 1 GeV, one can ignore inelasticities to a good approximation and, using unitarity, 
express in terms of phase shifts: 

N f^' , i X sin2(5;^(x) , 
W^{z) = — dxj ,,,, °^ a = 0,2 , 

' ^ ^ n ],Mi V X - 4M2 x{x - AMl){x - z) ' ^ 

In the upper part of the integration range, 0.5 GeV < ^/x < E, we can use the experimental 
phase shifts as discussed above. Below 0.5 GeV, experimental results for the phase shifts are 
rather limited. Fortunately, it is precisely in this range that we can make use of the chiral 
expansion of the amplitude. One expects that the chiral two-loop representation for the three 
phase shifts 6q, 6q and Si should be reasonably precise in the range 2M^ < y/s < 500 MeV: 
this expectation will be verified more quantitatively by comparing our results with those of a 
numerical solution of the Roy equations. We will therefore employ the chiral representation 
of the phase shifts in the dispersive integrals (|40|) below 0.5 GeV. When calculated in this 
way, the set of parameters depend on a, P and also, in a nonlinear way, on the four Aj's. 
A less precise procedure, but equally acceptable to the desired chiral order, is to use, below 
500 MeV, the imaginary parts of the S and P partial wave amplitudes in the chiral two-loop 
approximation (which is known to violate unitarity appreciably close to 500 MeV); i.e., one 
uses in the integrands of Eq. (|T5|) 



Im/o"(x) = J2wiix)lmK,ix) (a = 0,2) 

i 

lmfl{x) = {x- 4M^) ^ wi(x) ImKiix) . (41) 

i 

In this approximation, the dependence on Ai, A2 is simply linear and there is no dependence 
upon A3 and A4. In general, the system of equations ( p3D must be solved in a self-consistent 
way. 

Let us now discuss the numerical results. First, we consider the stability with respect to 
variation of the energy parameter E. Obviously, the parameters Ai, A4 should be independent 
of E. The validity of the formulas (|33|) , as discussed in the preceeding section, requires that E 
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should be smaller than, roughly, 1 GeV (see the discussion preceeding (p^) and at the same 
time that » 4M^. These conditions suggest to consider an interval 



500 MeV <E < 1000 MeV . 



(42) 



Variations of the results within this interval are shown in table 3 below ^ 



E (MeV) 


500 


600 


700 


800 


900 


lO^Ai 


-4.86 


-5.82 


-6.04 


-6.08 


-6.11 


lO^Aa 


9.68 


9.64 


9.60 


9.56 


9.55 


lO^As 


3.33 


2.55 


2.31 


2.20 


2.16 


10^A4 


-1.46 


-1.49 


-1.48 


-1.46 


-1.45 



Table 3: Results from the sum rules ( P^ j for several values of tie energy parameter E. The 

numbers correspond to a = 2 and (3 = 1.08 

Clearly, one does not observe exact stability. The variation is more significant for Ai and A3 than 
for the remaining two parameters, which are stable within 5%. A reasonable stability plateau 
forms for all four parameters for values of E above 700 — 800 MeV. Keeping in mind that 
E should remain sufficiently small such that the imaginary part of the D-wave contributions 
can be neglected, we shall assume in what follows that E = 800 MeV provides a reasonable 
compromise. 

Let us now discuss the results for the Aj's corresponding to various values of the two pa- 
rameters a and (3. This variation is of interest for the purpose of using the chiral formulas in a 
fit to determine a and (3 from experiment. Recall that a is expected to be close to 1 (within, 
say, 20%) according to standard chiral perturbation theory. In contrast, the generalized xPT 
can accomodate values of a as large as a = 4. At present, the best data available 
are compatible with a relatively wide range of values, 1 ^ a ^ 3. Some results, illustrating 
the dependence on a and /?, are collected in table 4 below (the numbers shown correspond to 
E = 800 MeV). 



a 


P 


lO^Ai 


10^X2 


10^ A3 


10^A4 


1.04 


1.08 


-5.68 


9.32 


2.21 


-1.48 


2.00 


1.00 


-6.60 


9.39 


2.12 


-1.33 


2.00 


1.08 


-6.08 


9.56 


2.20 


-1.46 


2.00 


1.13 


-5.70 


9.67 


2.26 


-1.54 


2.50 


1.08 


-6.27 


9.68 


2.19 


-1.45 


3.00 


1.08 


-6.45 


9.78 


2.19 


-1.44 


3.50 


1.08 


-6.62 


9.88 


2.18 


-1.44 


4.00 


1.08 


-6.79 


9.97 


2.17 


-1.43 



^As in I, we use F„ = 92.4 MeV and 



139.57 MeV. 
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Table 4: Results for the parameters Ai,...,A4 corresponding to several values of a and (3. 



The first line in the table corresponds to the values of a and f3 obtained in the standard xPT 
at chiral order 0{p'^) in the three fiavour case (see I). For the particular value a = 2, we have 
varied j3 in the range allowed by the Morgan-Shaw band. The table shows that as long as a 
and l3 remain in the range allowed by experiment, the variation of the values of the Aj's remains 
smaller than the error bars (see table 6). As a first approximation one may ignore this variation 
and adopt as average values 

lO^Ai ~ -6.1, 10^X2 ~ 9.6, 10^ A3 ~ 2.2, 10^ A4 ~ -1.45 . (43) 

These numbers update those given in our previous paper I 

lO^Ai ~ -5.3, 10% ~ 9.7, 10% ~ 2.9, 10% ~ -1.40 , (44) 

which were computed using a smaller value of E {E = 550 MeV), and with a slightly different 
treatment of the experimental data. The results are compatible within the uncertainties and 
the differences do not affect the numbers given in I for the threshold parameters. The results 
at higher energies are obviously more affected by small modifications of the values of Ai,...A4. 
For instance, computing the value of the phase of the e' parameter with the numbers from table 
4, one would find a value higher by 2 — 3° than the one given in I. 

Once the Aj's have been determined self-consistently, one can return to the sum rule ex- 
pressions and examine the contribution of each integration region. These are displayed in table 
5 which shows that the contribution from the resonance region, 0.5 GeV to 1 GeV, dominates 
the values of the three parameters Ai, A2 and A3, while the contribution from the low energy 
region is smaller though not negligible. The situation is different in the case of A4, for which 
the low energy region dominates. One may then fear that A4 could be overly sensitive to higher 
order chiral corrections. However, this does not seem to be the case: an estimate of the 0{p^) 
uncertainties is performed in the next subsection which shows that A4 is weakly affected (see 
table 6). 



(GeV) 


[2Af^,0.5] 


[0.5,1] 


[1,1.9] 


[1.9,00] 


103 Ai 


2.72 


-8.84 


0.14 


-0.10 


103 A2 


1.35 


7.69 


0.42 


0.10 


10^ A3 


0.49 


1.68 


0.03 


-0.001 


10^ A4 


-1.15 


-0.28 


-0.02 


-0.001 



Table 5: Contributions of four successive integration ranges to the values of Ai,...A4 with 

a = 2, /3 = 1.08. 
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3.3 Discussion of the errors 



We can identify two sources of uncertainties in the estimate of the Aj's: i) the errors af- 
fecting the experimental phase shifts and ii) the errors affecting the theoretical phase shifts, 
i.e. contributions of chiral order 0{p^) or higher that one might expect to become sizable in 
the neighborhood of 500 MeV. In order to estimate the uncertainties arising from the former 
source, we proceed as follows. First we take E to he E = 0.5 GeV. Using this value for E 
has the advantage that the experimental phases show up only in the integrals hi, h[ and not 
in the terms In computing we may use the approximation ( ^l]) so that the sum rules 
(|33|) can be solved analytically (since they reduce to a linear system) for Ai,...,A4. Finally, we 
assume that the errors on the phase shifts are piecewise constant as a function of energy (this 
assumption is seriously violated only in a narrow region around the KK threshold), and we 



have adopted the following values, an educated guess inspired by E^, |3^, |3T| : 



0.50 < < 0.95GeV 
0.95 < < 1.90 



4° 
9° 



A5\ 
A5\ 



1° 
1° 



A51 
A51 



T 
4° . 



(45) 



Under these assumptions the errors are given as simple analytic expressions in terms of a 
few phase shift integrals and it is a simple matter to calculate them. With the choice E = 
0.5 GeV, table 1 shows that most of the contributions to the integrals hi are concentrated 
in the integration region below 1 GeV. We have therefore neglected the errors coming from 
either the inelasticities or the partial waves with £ > 2, or from the asymptotic contributions. 
Finally, in order to estimate the error coming from the missing higher chiral orders in the low 
energy integrals, we have compared the result of calculating W^{z) using three evaluations of 
the integrands which have identical chiral expansions up to (and including) chiral order six, 
and which differ at higher chiral order: i) we use expressions ( ^OD together with the chiral 
approximation to the phase shifts (see formula (4.13) of I ) ii) we use the same expression 
but get the phase shift from a Fade approximant to the pertubative expansion which satisfies 
unitarity exactly (this generalizes to the 0{p^) case [] the idea used in the context of xPT in 
42| |) and iii) we use the chiral approximation (^iD , which violates unitarity. We have collected 



the results of these error estimates in table 6 where we show separately the contribution from 
each source, as discussed above. 





A5g 


A6l 


A6I 


0(/) 


Total 


lO^AAi 


1.02 


1.47 


0.40 


0.40 


2.23 


lO^AAa 


0.05 


0.26 


0.18 


0.20 


0.52 


lO^AAg 


0.34 


0.25 


0.14 


0.20 


0.64 


10^AA4 


0.01 


0.03 


0.09 


0.02 


0.12 



For the S and P partial wave amplitudes, the chiral expansion consists in three successive terms / 
/(2) + /(4) + f{6) and one finds fP^^'' = - /(4)//(2) + (/(4)//(2))2 _ /(6)//(2)]. 
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Table 6: Errors on Ai,...,A4 arising respectively from the experimental errors on 5q, 51 and 5q 
and from 0{p^) corrections to the amplitude below 0.5 GeV. 

The entries in the table were evaluated for a = 2, (3 = 1.08, but the errors show no 
significant variation with a or (3. We have added the errors associated with the experimental 
errors on the three phase shifts in quadrature, as it seems reasonable to assume that they are 
independent. The last piece of the error was added linearly. The reason why A2 and A4 have 
a small relative error compared to Ai and A3 is that the contribution from the isospin 1 = 
amplitude essentially drops out. (There is only an indirect small contribution which comes from 
the fact that the determination of A2 and A3 depends on Ai). The values of the parameters A2 
and A4 are thus to a large extent controlled by the p resonance. 

3.4 Comparison with former results 

The parameters A3, A4, being of chiral order 0{p^), have not been discussed previously in 
the literature. At O(p^), a potentially accurate determination of the parameters Ai and A2 
is possible, making use of the linear relationship with the low energy constants Li, L2 and 
L3 (see I: the relation is the same in the standard and in the generalized xPT). These three 
constants can be separately determined from the form factors of the Ku decay amplitude. On 
the theoretical side, this amplitude was analyzed at order one loop of the xPT ^ and 
estimates of higher order corrections were made |^5|. New high statistics data on decays 
would enable a fairly accurate determination of these constants. The best results available at 



present p5|, which make use of the experimental results of Rosselet et al. [O] together with 



the values of the tttt D-wave scattering lengths quoted in [3S], lead to 



Al = (-6.4±6.8)10"% A2 = (10.8 ±1.2) 10"^ (one loop) . (46) 



If one uses solely the constraints from the D-wave scattering lengths, then, according to |45 
the error on Ai would be increased by 50% and the error on A2 would be three times larger, 
while the central values would remain approximately the same. In comparing these central 
values with the numbers that we quote in table 3, one must keep in mind that in the latter 
Al and A2 include corrections of chiral order six. As far as the order of magnitude of these 
corrections is concerned, a reasonable guess should be provided by the values of A3 or A4. This 
leads one to expect that these corrections should not exceed 2 — 3%. On the contrary, we find 
that the 0{p^) corrections to the D-wave scattering lengths are significant. This is illustrated in 
table 7, where the values of a few threshold expansion parameters which pick up their leading 
contribution at O(p^) are shown, both at 0{p'^) and at O(p^) (the values of a, /3, Ai,...A4 being 
kept the same). 

These contributions explain why our central value for A2 is smaller than (^6[) by roughly 
10%. Apart from this effect, we believe that the errors on the LEC's (low energy coefficients) 
extracted from these D-wave scattering lengths have been somewhat overestimated by treating 
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them as independent experimental data. In reality, the numbers quoted in ||39| are obtained as 
sum rules from the Roy dispersive representation (|3^) projected on £ = 2, and are evaluated 
using experimental data at high energy and extrapolation of these data down to the threshold, 
based on the numerical solutions to the Roy equations of [^. It is obviously more efficient, as 
far as errors are concerned, to express directly the LEC's as sum rules. Concerning the P-wave, 
the results found in the chiral expansion for the threshold parameters a\, h\ and c\ tend to 
support the idea that the P-wave has a nearly pure Breit-Wigner shape down to the threshold. 
This indeed implies the following relations: 



4)^ 



(M^-4)' 



m-A) 



(47) 



where the mass and the width are expressed in units of the pion mass. One can check that 
these relations are rather well satisfied for a\ and h\ and remain qualitatively correct even for 
c\. An evaluation of h\ on the basis of sum rules was performed only recently as an outcome 



of new, rapidly convergent sum rules involving combinations of threshold parameters [^] and 
the result is 6| = (6 ± 4) x IQ-^. 





1-loop l-loop+2-loop Ref| 


39 




10^ al 

103 

Wh\ 


1.52 1.72 1.7 ±0.3 
0.20 0.14 0.13 ±0.3 
4.09 5.46 


W a\ 

lo^i 

10^ c\ 


0.30 0.58 0.6 ±0.2 
-4.79 -3.41 
-3.04 -3.54 
-1.82 4.64 



Table 7: Values of a few threshold parameters at one-loop and at two-loop accuracy (in 
units of Mt^+ ) for a = 2 and (3 = 1.08. The threshold parameters in the last four lines have no 

tree-level contributions at order O(p^). 

It is instructive to compare our results for the S and P partial-wave amplitudes derived in the 
two-loop approximation with a direct numerical solution of the Roy equations. This comparison 
allows to gauge the importance of 0{p^) terms which are present in the numerical solution and 
also allows to verify whether the driving terms agree. Numerical results for the phase shifts 
have been tabulated by Froggatt and Petersen corresponding to Oq = 0.30 and = —0.018. 



iQ — W.OW CllH^ U,Q 

Using the chiral expansion of the low-energy parameters (see appendix D of I) and the sum rules, 
we obtain the following numbers for the chiral parameters corresponding to these scattering 
lengths: a = 2.84, /? = 1.09, lO^Ai = -6.32, lO^Aa = 9.77, 10^ A3 = 2.20 and 10^A4 = -1.46. 
Figure 1 shows that the phase shifts computed from the two-loop chiral expansion agree rather 
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well with the results of up to an energy ~ 500 MeV, thereby justifying the assumption 



made in the derivation of the sum rules to trust the chiral representation precisely in this range. 

3.5 Implications for the standard two-loop amplitude parameters 

While the present article was being completed, the computation of the vrvr scattering am- 
plitude to two-loop accuracy has been achieved in the framework of the standard xPT in 
Ref. [^. Comparing the expression of A{s\t,u) obtained by the authors of fl^ with the one 
we had derived in I, we conclude that they coincide up to order 0{p^) contributions, provided 
one identifies the six constants bi which appear in the amplitude of [jl9| with our parameters 
a, P, Aj as follows |: 

1, Ml f^^ 1 103 

^3 = Ai + -A2 - 6A3 + 



2 F3 \ 122887r2 46087r^ 



1, Ml 
:A2 + 



^ 2 " 15367r4F2 

1 19 
^5 = A3 — -A4 



4 15367r4 

- ~ 4608^ 

bi = - 1) - oT7^(/5 - 1) + 4Ai - 8A3 + 



?>My ' " FIK" 11527r2 1447r4 

/ 19 1 \ 

62 = ±ZL(/3_i)_4A, + l!:i^ I2A3 + — t)- (48) 
MY Fl\ 184327r2 967rV ^ ' 

From the last two formulae, it is apparent that in the standard framework, both a and have 
to become equal to one in the chiral limit. (In G^PT, a stays away from 1 and b\ becomes 
large, signaling the breakdown of the standard expansion.) 

From our preceding determination of the four constants Aj, we can now deduce the values 
of the four constants 63, 64, 65 and b^. In order to make the result as accurate as possible, the 
0(j9^) values of a and /? in the standard framework would be required. Unfortunately, we see 
no way to extract these values from the content of [1^. However, we expect that a and /5 will 



remain close to their leading order values a = /? = 1 in the standard framework. The Oijf) 
values (in the two flavour standard xPT) a = 1.06, (3 = 1.10 provide already some information 
about the expected corrections. Using these values, we obtain (taking = 93.2 MeV here in 



order to conform with [p!9| ]): 



63 = (-3.7 ± 2.4) 10"^ 64 = (4.8 ± 0.3) 10"^ 

65 = (1.4 ±0.6) 10"^ 66 = (1-OiO.l) 10^^ (49) 

If one assumes that the 0{p^) corrections to a and P are of the order of the square of the O(p^) 
corrections, then the above numbers are not modified. 



^We are indebted to G. Colangelo for pointing out a mistake in the first version of the manuscript. 
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4. Summary and conclusions 

i) Among the six parameters a, /3, Ai,...,A4 that define the two-loop tttt scattering am- 
phtude, the parameter a is the most sensitive to the chiral structure of the QCD vacuum. 
Furthermore, it is relatively weakly affected by higher orders of xPT, less affected than, for in- 
stance, the S-wave scattering lengths. For these reasons, the parameter a represents a suitable 
quantitative characteristic of the strength of quark condensation, and it should be determined 
experimentally in order to disentangle the large and small condensate alternatives. To accom- 
plish this, one needs new high-precision low-energy vrvr scattering data and a model-independent 
determination of the parameters Ai, A4. Indeed, this would allow a measurement of a (and 
j3) from a simultaneous fit to several low-energy observables such as 5q — 5\ |^ and a^ — a^ ||21|| . 



The present paper contributes to this program by providing a rather accurate determination of 
the parameters Ai, A4 using existing tttt scattering data at medium energies. 

ii) The method of determination follows from a rather systematic procedure: first, one 
establishes a dispersive representation for the low-energy amplitude in terms of absorptive 
parts and two unknown subtraction constants. This representation is displayed in Eg . (^0|) , 
and it holds up to corrections of order 0[(p/i?)®], with 4M^ << E'^ 1 GeV^. Its derivation 
parallels the derivation of the Roy equations, using twice-subtracted fixed-t dispersion relations 
and exploiting crossing symmetry. The second step consists in identifying, up to and including 
0(jP) accuracy, the dispersive representation (|30| ) and the perturbative two- loop amplitude in 
a whole low-energy domain of the Mandelstam plane. The comparison yields six, and only six, 
sum rules for a, /3, Ai, A2, A3, A4. The sum rules for a and [3 are only barely convergent 
and they were ignored. One remains with four rapidly convergent sum rules for Ai, A4. Our 
method minimizes the error in the determination of these parameters, since it makes full use 
of the information contained in crossing symmetry. 

iii) In evaluating the sum rules (0), one uses the existing vrvr data in the energy range 
0.5 GeV < ^/s < 1.9 GeV which were extracted from unpolarized high statistics ttN — >■ 
TTT^N high-energy production experiments. The outcome is rather sensitive to these data; in 
particular, their error bars constitute the main source of uncertainty in the values of the Aj's (see 
table 6). Our machinery is ready to accept as input any other set of medium energy vrvr phases 
and inelasticities, provided they are consistent with the Roy-type dispersion relations. (The 
standard vrvr phases extracted from the old Cern- Munich experiment have been recently 



criticized ^^l- However, no alternative phases have been proposed.) The high-energy tails in 
the sum rules ( |33D are estimated using the Regge-pole model and they are found to be negligible. 
Finally, the contribution of the low-energy range 2Mn- < y/s < 0.5 GeV is evaluated using the 
perturbative two-loop amplitude itself. This contribution introduces a weak dependence of 
the parameters Ai,...A4 on a and /? illustrated in table 4. This low-energy part of the sum 
rules encompasses the bulk of infrared contributions represented by the chiral logarithms in an 
explicit xPT calculation. 

iv) Our final result can be read off from tables 3, 4 and 6. The variation of the A^'s with the 



22 



cutoff E exhibits a plateau for E ~ 700 — 900 MeV. Choosing for the central values E = 800 
MeV as well as the center of the variations with respect to a and to (3, we obtain as the final 
result of this paper 

lO^^Ai = -6.1 ±2.2 , 10^ As = 9.6 ±0.5 

10% = 2.2 ±0.6, 10% = -1.45 ±0.12 . (50) 

The central values may be further refined taking into account the small dependence on a and 
on (3. In any case, the variation with a and (3 is smaller than the error bars displayed in (|50[). 
Giving to the parameters Ai,...,A4 the values (|50|) , the two-loop perturbative amplitude becomes 
a faithful analytic low-energy representation of the numerical solution of the Roy equations, 
which were used to establish the values, usually quoted as "experimental", of tttt threshold 
parameters [^]. Indeed, for each pair of values of a and (3 ( Oq and ) the perturbative 
amplitude reproduces the S-wave slopes, and the P, D and F-wave threshold parameters as 
quoted in ||39 . 



v) The parameters Ai,...,A4 are determined globally, without any reference to their ex- 
pansion in powers of quark masses and chiral logarithms. Consequently, the resulting values 
(|50| ) concern both the standard and the generalized xPT (modulo a small variation with a 
and f3 displayed in table 4). It would be interesting to see, whether the "chiral anatomy" of 
these parameters resulting from the recent standard xPT two-loop calculation can be used 
to predict the values of Ai,...,A4 (or, equivalently, 63,...,6e) in agreement with their present 
determination. 
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Figure l:Comparison of the S- and P- wave phase shifts in the chiral expansion at two- loop 
with the results of a numerical solution of the Roy equations quoted by Froggatt and Petersen 
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